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I. INTRODUCTION 



A. TAYLOR- GORTLER VORTICES 

Since the early part of the twentieth century, consider- 
able research has shown that fully developed laminar flow 
along a concave wall does not remain two-dimensional [Refs. 

1, 2, 3]. The flow instead forms a system of spiral vortices, 
of counter rotating pairs, whose axes are in the principle 
flow direction. This phenomenon is the result of the varia- 
tions in the centrifugal forces acting on the fluid particles, 
and is known as Taylor-Gortler vortices. Figures 1 and 2 
illustrate the type of fluid motion just described. 

In a channel that is curved in the streamwise direction, 
those fluid particles located in the center of the flow cross- 
section are subjected to higher centrifugal forces than those 
fluid particles traveling along the channel's boundary wall. 

As a result, the fluid in the center of the channel moves out- 
wardly towards the concave boundary. As the process continues, 
the fluid particles near the boundary wall, move in a spanwise 
direction, and finally radially inward replacing the outwardly 
moving particles. These particles then come under the same 
centrifugal forces and the process repeats itself. The re- 
sulting cyclic motion, causes the formation, and propagation, 
of the counter-rotating Taylor-Gortler vortices which are con- 
sidered primarily a laminar flow phenomenon. 
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Figure 1. Schematic of Taylor vortices between cylinders. 




Figure 2. 



Schematic of Taylor - Gortler vortices 
channel . 



in a curved 
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It has been observed [Ref. 4] that the heat transfer rate 
from flow along a concave curved wall is greater than that 
for flow along a straight wall. It is thought that the addi- 
tional mixing provided by the secondary motion of the Taylor- 
Gortler vortices may account for this. It is also thought that 
the cross-hatching that has been observed in reentry vehicles 
can be explained in part by the presence of streamwise vortices 
similar to the Taylor-Gortler vortice phenomenon. 

There are many possible applications that could result from 
a more thorough understanding of Taylor-Gortler vortices and 
their effect on heat transfer and fluid flow characteristics. 

By taking advantage of the increased heat transfer rate, im- 
proved heat exchanger designs and improved turbine blade cool- 
ing could result [Refs. 5, 6, 7]. 

B. HISTORY 

The instability of an inviscid fluid flowing past a curved 
boundary, was first considered by Lord Rayleigh in 1916 [Ref. 

8], By assuming that the fluid was non-viscous, he determined 
that for the motion of an inviscid fluid to remain stable, its 
circulation must increase with increasing radius. G. I. Taylor, 
in 1923, [Refs. 1, 9], continued this study with an extensive 
analytical and experimental study of viscous fluids. His in- 
vestigations focused on the flow between two cylinders. For 
the case where the inner cylinder was rotated while the outer 
cylinder was held stationary, Taylor ascertained that such 
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flows become unstable when the value of a dimensionless para- 
meter exceeded a critical value of 41.3. The parameter, 
known as the Taylor number is defined as: 

Ta = Re /jC 



where T d ' is the width of the gap, assumed small when compared 
to ' Ri ' , the radius of the inner cylinder, and 'Re' is the 
Reynolds number. The Reynolds number is based on the peri- 
pheral velocity of the inner cylinder. Taylor determined 
that for those cases in which the value of the Taylor number 
exceeded 41.3, a secondary motion developed and the Taylor 
vortices formed. 

Instability of a similar nature is also observed when a 
viscous fluid flows in a curved channel due to a pressure 
gradient acting along the channel wall. This problem was first 
considered analytically by W. R. Dean [Ref. 10] in 1928, for a 
channel formed by two concentric cylinders, where the radius of 
the innder cylinder was large in comparison to the spacing be- 
tween the inner and outer cylinder walls. Dean concluded that 
there would be an initiation of the flow instability, and the 
propagation of vortices, when a similar dimensionless parameter, 
the Dean number, exceeded a value of 36. The Dean number is 
defined as: 



De 



Re 
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where 'd' represents the channel half-width, ' Ri ' is the inner 
cylinder radius, and 'Re' is the Reynolds number based on the 
mean velocity of the undisturbed flow. The analytical work 
of Dean was later verified by W. H. Reid [Ref. 11], using an 
approximate solution. 

In 1940, H. Gortler studied the stability of laminar bound- 
ary layer profiles on curved walls under the influence of small 
disturbances. He found that these disturbances were similar 
to the vortices studied by G. I. Taylor. By approximate numer- 
ical calculations, Gortler concluded that the disturbances or 
vortices were produced only on the concave boundary walls. He 
also concluded that the overall flow profile appeared to remain 
laminar in nature [Ref. 2]. These results were verified, with 
an exact solution, by G. Hammerlin, and reported by H. Schlichting 
[Ref. 12] in 1955. A. M. 0. Smith completed an even more ex- 
tensive numerical analysis that further substantiated these 
findings [Ref. 3]. Recently with the use of hot wire anemometry, 
laser doppler systems, and flow visualization techniques, [Ref. 
13, 14], the result of these numerical solutions have been 
demonstrated experimentally. In 1976, Y. Aihara [Ref. 15] con- 
ducted a non-linear analysis of Gortler vortices. 

As interest began to develop concerning the effects of the 
secondary flows associated with the Taylor-Gortler vortices, 
many studies were published concerning the influences of these 
vortices on the transfer of heat in the laminar and turbulent 
flow regimes. In 1955, F. Kreith [Ref. 16], studied the 
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influence of heat transfer with respect to the curvature of 
the boundary wall for fully turbulent flows. He concluded 
that the heat transfer from the heated concave boundary wall 
was considerably higher than that transferred from the convex 
boundary wall of the same curvature and under similar turbu- 
lent conditions. 

In 1965, L. Persen [Ref. 17], considered the special cases 
of very high and very low Prandtl number fluids. He related 
the increase in heat transfer from a curved wall to the pre- 
sence of the Taylor-Gortler vortices. 

There has been only a limited amount of published litera- 
ture dealing with the flow and heat transfer in curved channels 
of rectangular cross-section and large aspect ratios. Much 
of what has been published involves the development of numer- 
ical approximations and solutions for heat and mass transfer 
in curved ducts of various geometries. K. Cheng and M. Akiyama 
[Ref. 18] developed a numerical solution for forced convection 
heat transfer with laminar flows in curved channels of rec- 
tangular cross-section. However, they were concerned with 
cross-sections with small aspect ratios. 

In 1976, A. A. Shibani and M. N. Ozisik [Ref. 19] solved 
the heat transfer problem between parallel plates with tur- 
bulent flow, for the case of uniform wall temperature. They 
used matched asymptotic expansion techniques for a wide range 
of Prandtl numbers. 
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Other experimental and analytical studies that are worth 
citing with regard to this present study follow. Y. Mori [Ref. 
20], obtained results for hydrodynamical ly fully developed 
flows with constant wall heat flux, in curved channels of 
square cross-section. W. M. Kays and E. Y. Leung [Ref. 21] 
reported solutions for turbulent flow heat transfer in a con- 
centric circular tube annulus with a fully developed velocity 
profile and constant heat rate per unit length, for a fluid 
of Prandtl number 0.7. Results for large aspect ratio channels 
with rectangular cross-sections for laminar flows were reported 
by M. Durao [Ref. 22] and J. Ballard [Ref. 23], while R. Holihan, 
Jr. reported results for laminar and transition flows [Ref. 

24]. P. F. Brinich and R. W. Graham [Ref. 25], reported 
results for turbulent flows in a rectangular curved channel 
with an aspect ratio of 6, for the inner wall heated, the 
outer wall heated, and both walls heated. 
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II. INTENT OF THE STUDY 



The purpose of this investigation was to examine the 
effect of streamwise curvature on the heat transfer rate in 
a curved rectangular duct of large aspect ratio. The flow 
velocities examined were in the turbulent regime. The results 
of this study were compared to the heat transfer rates for 
the same range of velocities, in a straight duct of identical 
aspect ratio. 

Taylor-Gortler vortices are considered to be a laminar 
flow phenomenon. Enhanced heat transfer in curved channels 
for transition and turbulent flows has been observed. Even 
though it has not been proven that the vortices continue to 
propagate at the higher flow velocities, this enhancement in 
the transfer of heat has been attributed to the secondary 
flow velocity components of the Taylor-Gortler vortices. It 
is believed that these secondary components transport the 
heated fluid from the outer concave wall, inward toward the 
opposite wall of the channel, displacing the cooler fluid 
particles and causing them to move toward the heated concave 
wall. It was expected that similar results would be observed 
in this study. 

This investigation was conducted using a single channel 
with a rectangular cross-section and constant aspect ratio. 

The channel incorporated both a straight test section and a 
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curved test section. The results obtained in each test section 
at approximately the same flow rates were compared in an effort 
to determine the effects of the Taylor-Gortler vortices on the 
transfer of heat. Also, the results of both the straight and 
curved sections were compared to the results of Ballard [Ref. 
23] and Holihan [Ref. 24] . 

The straight section results of this study were compared 
to the results of Shibani and Ozisik [Ref. 19] for turbulent 
flow between parallel plates. In addition, the results of 
Kays and Leung [Ref. 21] and the Dittus -Boelter equation 
[Ref. 26] were compared for Reynolds numbers in the range of 
this study. Kays and Leung's data was for turbulent flow in 
annular passages while the Dittus -Boelter equation was for 
turbulent flows in circular tubes. 

The curved section results of this study were compared 
with the results of Brinich and Graham [Ref. 25] for turbulent 
flow in a rectangular curved channel. 
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III. EXPERIMENTAL WORK 



A. DESCRIPTION OF THE APPARATUS 

To meet the objectives of this investigation, a channel 
of rectangular cross-section was used. Details of its con- 
struction are described in references 22 and 23. Plexiglas, 
0.635 centimeters thick, formed the inner and outer walls of 
the channel. These walls were separated by 0.635 centimeter 
spacers that also served as the sides of the channel. 

The channel, shown in Figure 3, had a straight section, 
122.0 centimeters in length, followed by a curved section of 
180 degrees of arc, followed by another straight section, 91.4 
centimeters in length. The inner curved wall of the channel 
had a radius of curvature of 30.5 centimeters. The rectan- 
gular channel was .635 centimeters high and 25.4 centimeters 
wide resulting in an aspect ration of 40. The hydraulic dia- 
meter was 1.229 centimeters. The cross-sectional area of the 
channel was 16.13 square centimeters. The wetted perimeter 
was 52.07 centimeters. Figure 4 shows a cross-sectional view 
of both the straight and curved sections of the channel. 

An entrance bell constructed of plexiglas was connected 
to the straight section of the channel. It was designed and 
manufactured according to ASME nozzle standards, with an 
elliptical curved base on a major axis equal to ten inches 
and a minor axis of one inch. Cheese-cloth was attached to 
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Figure 3. Illustration of equipment and apparatus. 
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Figure 4. Cross-sectional view of channel. 



the entrance nozzle to prevent foreign matter from entering 
the flow channel. An aluminum exhaust nozzle was attached 
to the exit of the channel and directed the flow from the 
channel into tivo inch diameter pvc piping. The tx^o inch 
pvc piping contained a concentric orifice with a diameter 
of 1.065 inches and was constructed according to the ASME 
Power Test Code. Pressure taps on either side of the orifice, 
at 1 and 1/2 pipe diameters from the orifice, and connected 
to manometers were used to calculate the mass flow rate, of 
the working fluid, which was air at room temperature. 

The flow of air was drawn through the channel by an 
electrically driven Spencer Turbo Compressor, 30 horsepower 
at 3500 rpm, rated at 550 cubic feet per minute at 70 degrees 
fahrenheit and one atmosphere. The channel contained two 
test sections from which experimental heat transfer data 
could be obtained. The straight test section was located at 
a distance downstream of the entrance bell that would ensure 
hydrodynamically fully developed flow. The straight test 
section was 29.2 centimeters in length with a heated test 
section area of 741.7 square centimeters. The curved test 
section was located in the lower half of the curved portion 
of the channel. It was 28.3 centimeters in length, subtend- 
ing an arc of 53.1 degrees, and having a heated test section 
area of 718.4 square centimeters. 

In each of the heated test sections, Temsheet, a carbon 
impregnated porous paper with a uniform electrical resistivity. 
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was glued to the outer wall of each section. Joulean heating 
was used to heat the flow of air through the channel. Since 
the electrical resistance of the Temsheet is not constant, 
but varies slowly with temperature, a precision resistor with 
an electrical resistance of 2.0262 ohms was connected in series 
with the Temsheet to allow the calculation of the instantaneous 
power being supplied. Detailed schematics of each test section 
are shown in Figures 5 and 6. 

The variables that were measured and used in this invest- 
igation were: 

(1) the temperature of the air entering the channel (T^ ) 

(2) the temperature of the air leaving the channel (T Qut ) 

(3) the temperature of the heated outer wall for each 

test section (T WQ ) 

(4) the temperature of the unheated inner wall for each 
test section (T .) 

(5) the temperature between the three layers of insulation 
at each test section (T^ ng ) 

(6) the voltage across the precision resistor (V ) 

(7) the voltage across the Temsheet or heater (V^) 

(8) the temperature of the air at the orifice C T q ^ ) 

(9) the pressure upstream of the orifice (P ) 

Lip 

(10) the difference in the upstram and the downstream 
pressure across the orifice (AP) 

(11) the atmospheric pressure (P atm ) 
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Figure 5. Straight test section, detailed schematic. 
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Figure 6. Curved test section, detailed schematic 
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Each of the temperature measurements was made with Copper- 
Constantan thermocouples. Those thermocouples that could be 
removed easily from the channel were calibrated using a 
ROSEMONT Commutating Bridge model 920A, and a ROSEMONT Constant 
Temperature Bath model 913A. A second order polynomial was 
fitted to each set of data. The resulting coefficients dif- 
fered by less than 1.0 percent, so an average of each coeffi- 
cient was used in converting emf to degrees centigrade in the 
reduction of the data obtained. 

A total of sixty-one Copper-Cons tantan , 30 gauge thermo- 
couples were located throughout the test apparatus to record 
the desired temperatures. Four thermocouples were connected 
in parallel to read an average temperature at the entrance 
and exit to the channel. Five thermocouples were connected 
in parallel and inserted between each of the layers of insula- 
tion at each of the test sections. In addition, individual 
thermocouples were positioned in each of the test sections, 
eight in direct contact with the Temsheet and eight in the 
plexiglas of the unheated inner wall. Small diameter holes 
in the plexiglas were drilled to allow the installation of 
the thermocouples. One thermocouple was inserted in the 
two inch pvc pipe downstream of the orifice in accordance 
with the ASME Power Test Code to record the temperature of 
the air flowing through the orifice. The beads of the thermo- 
couples in contact with the Temsheet were electrically insu- 
lated from the Temsheet with ENMAR Heat Resisting Glyceryl 
Phthalate . 
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The thermal insulation consisted of 1/4 inch layers of 
Armstrong Armaflex 22 Sheet Insulation, a flexible foamed 
plastic material, on the heated side of each test section, 
and 1/2 inch Armaflex 22 insulation elsewhere on the test 
apparatus. Three layers of the 1/4 inch insulation were 
used and covered an area slightly larger than each of the 
heated areas. By positioning thermocouples between these 
insulation layers, the heat lost to the environment could 
be computed at each test section. The heat loss through 
the 1/2 inch insulation surrounding the rest of the channel 
due to the heated air flowing through the channel was small 
in comparison. See Figure 7 for a detailed sketch of the 
thermocouple placement and attachment. The insulation was 
held in place by the use of Adhesive Heat Resistant ven- 
tilation Duct Tape. The entire channel and all connections 
between the channel, the pvc piping, and the orifice were 
sealed with General Electrical Silicone Rubber Sealant Caulk 
to ensure there was no leakage of air into the channel or the 
piping which could affect the temperatures recorded as well 
as the mass flow rate. Two aluminum electrodes, 0.318 centi- 
meters thick, were inserted in between the Temsheet and the 
plexiglas in each section as shown in Figure 7. 

The data acquisition system used for this investigation 
was a Hewlett Packard 3054A Automatic Data Acquisition/Control 
System consisting of a 3456A Digital Voltmeter and a 3497A 
Data Acquisition/Control Unit. Also used in conjunction with 
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Figure 7. Thermocouple placement in test sections. 
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Figure 8. Photograph of test channel. 
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Figure 9. Photograph of flow measuring apparatus. 
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Figure 10. Photograph of data acquisition system. 



this data acquisition system, were a Hewlett Packard 9326 
Computer terminal and 2671G Printer. Pressure measurements 
were taken from Meriam vertical manometers, one whose fluid 
was water with a 0 to 60 inch range, and the other whose 
fluid was mercury, calibrated to read inches of water with 
a range of 0 to 415 inches. A photograph of the channel and 
associated test equipment are shown in Figures 8, 9, and 10. 

B. EXPERIMENTAL PROCEDURES 

Experiments were conducted by varying the mass flow rate 
of air corresponding to Reynolds numbers from approximately 
12000 to 23000, where the Reynolds numbers are based on hy- 
draulic diameter. For each mass flow rate, data was taken 
automatically by the data acquisition system and the data 
was immediately reduced and printed for examination and com- 
parison. The experimental procedures followed were the same 
for the straight and curved test sections. The results ob- 
tained from the straight test section served as a baseline 
for the comparison of the curved section results. Prelimi- 
nary runs were performed to determine the time required for 
the test rig to come to steady state. For this experiment 
three hours were allowed for steady state to be reached. 
After three hours, data was taken at ten minute intervals 
to observe changes in parameters to ensure that a steady 
state condition had been achieved. The criteria for steady 
state was based on three variables: 
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(1) the mass flow rate 

(2) the heated boundary wall temperature 

(3) the heater voltage 

When these variables varied by less than two percent over a 
ten minute interval, it was considered that a steady state 
condition had been achieved. 

It was also necessary to determine the approximate power 
setting required to bring the heated wall of each test section 
to approximately 50 degrees Centigrade. It had been deter- 
mined by Holihan [Ref. 24] that this temperature was sufficient 
to ensure a 20 degree Centigrade difference in the heated \vall 
and the unheated wall. 

For each run, the atmospheric pressure, the pressure dif- 
ference across the orifice, and the pressure upstream of the 
orifice were entered into the computer for the mass flow rate 
calculations. All other data was acquired by the data acquisi- 
tion system directly from the experimental apparatus. 

To compute the instantaneous power supplied the following 
relationship was used: 



Qp 



V H V 
R PR 



R 



The precision resistor voltage and the heater voltage were 
read by the data acquisition system and the resistance of the 
precision resistor was a known constant. A diagram of this 
circuit is shown in Figure 11. 
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Figure 11. Circuit used to measure power supplied to 
the Temsheet. 
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IV. PRESENTATION OF DATA 



A. ANALYSIS 

The constant heat flux surface for the heated wall of each 
test section was approximated by the uniform electrical resis- 
tivity of the Temsheet. The insulated unheated wall was con- 
sidered adiabatic since the heat losses through that wall were 
negligible. The channel’s large aspect ratio of 40 allowed 
the modeling of the channel as parallel plates. The straight 
portion of the channel upstream of the straight test section 
was of sufficient length to ensure that the flow was hydro- 
dynamically developed for the flow velocities of this study. 
The straight portion of the channel downstream of the curved 
test section was also of sufficient length to ensure that the 
flow exiting the channel was thoroughly mixed and that the 
temperature, T QUt , was an average bulk temperature. 

Based on these assumptions, the experimental configura- 
tion was modeled as forced convection between parallel plates, 
for hydrodynamically developed and thermally developing flows, 
subjected to a constant heat flux. The boundary conditions 
are one wall at a constant heat flux and the other wall 
insulated . 

To analyze this problem, several quantities were defined 
as follows: 

The heat convected to the air was calculated using the 
express ion : 
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^air ™ ^p ^out 



where 'C - ' was the specific heat of the air at constant 

pair r 

pressure and 'm' was the mass flow rate of the air. 

The average heat transfer coefficient between the heated 
wall and the flow of air in the channel was defined by the 
equation: 



h = 



^air 

a pl at 



where 'Q a ^ r ' is defined above, 'Ap^* was the area of the 
Temsheet in the test section, and ’AT' was the difference 
between the average heated wall temperature (T ) and the 
average bulk temperature (T^^) • The average bulk tempera- 
ture was defined as the arithmetic mean of the entrance and 
exit temperatures (T in> T Qut ) . 

The average Nusselt number was then calculated as follows 

»- * D hd 

Nu = 

K • 
air 



In this equation ' D^’ is the hydraulic diameter and ' K a ^ r * 
is the thermal conductivity of air. 

The Reynolds number was calculated for each test run as 
follows : 



m D 



Re 



hd 



hd 



A c ^air 
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where again 'm' and 'D^' are the mass flow rate and hydraulic 
diameter of the channel, 'A ' is the cross-section of the 
channel, and 'y a ^ r ' is the dynamic viscosity of the air. 

For the curved section runs a Dean number was defined as: 




where 'Re' is the Reynolds number based on hydraulic diameter, 
'Dj^' is the hydraulic diameter, the 'Ri' is the radius of 
curvature of the unheated inner convex wall. 

A sketch of the control volume and a set of sample calcu- 
lations for one test run of the curved section are given in 
Appendix B. 

B . RESULTS 

The data obtained from each experimental run was evaluated 
utilizing the expressions described in Part A above. The major 
parameters resulting from this evaluation are shown in Tables 
I and II. Table I contains the straight section results and 
Table II contains the results of the curved section. A plot 
of the average Nusselt number verses Reynolds number is given 
in Figure 12 for the comparison of the results. Error bands 
have been indicated as a result of an uncertainty analysis. 

A sample calculation for the uncertainty analysis is given in 
Appendix A. 

The results indicate an increase in the rate of heat 
transfer with increasing Reynolds number for both the 
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Figure 12. Straight vs. curved section results, present study. 



TABLE I 



SUMMARY OF STRAIGHT TEST SECTION RESULTS 



Re hd 


Q • 

x air 


E 


AT 


Nu 




(W) 


(W/m C) 


(C) 




12940 


114.43 


63.55 


24.27 


29.9 


14520 


113.55 


67.05 


22.82 


31 . 5 


17020 


120.00 


76.70 


21.09 


36.1 


20480 


117.37 


84.30 


18.76 


39.7 


22980 


121 . 31 


92.04 


17.76 


43.4 


Reynolds 


number and 


Nusselt number 


are based 


on hydrau 



diameter . 
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TABLE II 



SUMMARY OF CURVED TEST SECTION RESULTS 



Re hd 


De 


Q • 

H air 

(W) 


13490 


2707 


115.28 


16950 


34 0 3 


117.50 


21140 


4243 


114 . 94 


23320 


4680 


113.43 



E 


AT 


Nu 


(W/m C) 


CCD 




77.23 


20.79 


36.2 


91.26 


17.93 


42.9 


105.46 


15.18 


49.7 


111.47 


14.17 


53.9 



Reynolds number, Dean number, and Nusselt number are 
based on hydraulic diameter. 
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straight and curved test sections. In addition, the heat 
transfer rate was higher for each Reynolds number investigated 
in the curved section. 

In the investigation by Holihan [Ref. 24] a straight line 
least squares correlation of his data resulted in: 



Nu = 0.983 Re d °' 25 
for laminar flows, and 


300 < Re d < 2000 


Nu = 0.110 Re^’ 8 ^ 


2000 < Re, < 3000 
a 



for transition flows in the straight test section. The curved 



section data yielded: 




Nu = 0.065 Re d ^’^ 7 
for laminar flow, and 


300 < Re d < 2000 


Nu = 0.117 Re d 0,58 


2000 < Re d < 3000 


for transition flows. 





The data for the present study resulted in 



Nu = 0.063 Re hd °’ 65 


10 4 < Re^ d < 2.5 x 10 4 



for turbulent flow in the straight test section, and 
Nu = 0.040 Re hd 0 ' 72 10 4 < Re hd < 2.5 x 10 4 

for similar flow in the curved section. 
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In the turbulent regime of this study, the rate of heat 
transfer was about twenty percent higher in the curved section 
when compared with the straight section. This compares favor- 
ably with earlier studies by Kreith [Ref. 16] , where he re- 
ported an increase in the rate of heat transfer along a 
concave wall of from twenty-five to sixty percent for 

Reynolds numbers, based on hydraulic diameter, between 
4 6 

10 and 10 . Holihan reported an increase of fifteen per- 
cent for laminar flow and thirty percent for transition 
flow [Ref. 24] . Ballard reported an increase of eleven 
percent in the heat transfer rate for laminar flows [Ref. 23]. 
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V. DISCUSSION AND CONCLUSIONS 



Holihan [Ref. 24], in his study, determined that there 
was a negligible difference in the fluid bulk temperature 
and the unheated wall temperature. He additionally demon- 
strated that the radiated heat transfer from the unheated 
wall surface was minimal. These results were verified for 
the present study and provide the basis for the assumption 
that the heat transfer to the air flowing through the channel 
was solely by convection from the heated Temsheet. 

As was mentioned earlier, the high aspect ratio of the 
channel provided the basis for the assumption that the ex- 
perimental apparatus, as configured, could be modeled as 
infinite parallel plates. Holihan's data [Ref. 24] also 
tended to verify this assumption. His experimental data, 
in the laminar flow region, approached the theoretical limit 
for average Nusselt number for parallel plates with one wall 
heated at a constant heat flux and the opposite wall 
adiabatic [Ref. 27]. 

Based on the assumptions mentioned above, comparisons were 
made with other analytical and experimental results that were 
of the same or similar problem. Again the problem being, flow 
between infinite parallel plates with one wall at a constant 
heat flux and the opposite wall adiabatic. 
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A comparison of the experimental results of this study 
and the experimental results of Ballard [Ref. 23] and Holihan 
[Ref. 24] are shown in Figure 13, for the straight section, 
and Figure 14 for the curved section. Figure 15 shows a com- 
pilation of all the data on one plot. Ballard's data was 
strictly in the laminar region, while Holihan' s data covered 
the laminar as well as transition regions. The channel aspect 
ratio was 40 for each of these studies and the experimental 
procedures were similar in each case. Ballard and Holihan 
compared their data to the analytical studies by McCuen, 

[Ref. 28] , for heat transfer between infinite parallel plates 
with constant wall temperatures and heat flux. Their results 
plotted above the analytical solution and the difference was 
attributed to the difference in the geometries used in the 
studies, and to the limitations inherent in any experimental 
work. The side wall effects and the inability to totally 
account for all the heat transfer processes and/or losses 
were also factors that may have caused the difference. 

A comparison of analytical and experimental results with 
those of the present study for the straight test section is 
shown graphically in Figure 16. The results of Shibani and 
Ozisik [Ref. 19] is shown for turbulent flow between parallel 
plates. The correlation used was 

Nu = 12 + .024 Re * 824 
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4 6 

for 0.1 < Pr < 10 and 10 < Re < 10 . The Dittus -Boelter 

equation [Ref. 26] for heat transfer in a straight tube for 
constant wall temperature is plotted from the following 
correlation : 

Nu = .02 Re* 8 

for a Prandtl number equal to 0.71. Additionally, Kays and 
Leung [Ref. 21] obtained data for heat transfer in annular 
passages, with one wall heated at a constant heat rate. 

The data shown is for r* = 1.0 which equates to parallel 
plates . 

The curve section data is plotted in Figure 17 and is 
compared with the experimental results of Brinich and Graham 
[Ref. 25] for flow and heat transfer in a curved channel. The 
accuracy of the data points used are subject to errors, in 
that the actual values were not given in their study and had 
to be taken from a plot of Stanton number versus Reynolds 
number as given in the reference. 

The data from the present study fell below the correlations 
that were used as comparisons. These differences can also be 
attributed to the differences in geometries and conditions of 
heating used. Experimental error and side wall effects may 
have also contributed to the overall difference in results. 

The data from the present study and that from the Brinich 
and Graham study seem to compare favorably and exhibit a simi- 
lar trend. Irrespective of the differences noted above, it 
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Figure 17. Comparison of present data with Brinich 
and Graham, for turbulent flow, curved 
section . 



appears that there is a definite enhancement of the heat 
transfer rate in the curved test section for laminar, tran 
sition, and turbulent flows. Streamwise curvature, result 
ing in secondary flows, seems to cause the enhancement. 
Whether the Taylor- Gortler vortices continue into the 
transition and turbulent flow regimes remains to be seen. 
Continued study of the effect of streamwise curvature and 
the development and propagation of the Taylor-Gortler vor- 
tices throughout all flow regimes is required. 
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VI. RECOMMENDATIONS 



To better understand the effects of streamwise curvature 
on the enhancement of heat transfer, a significant amount of 
study still remains. With the introduction of today's data 
acquisition systems, in conjunction with portable computers, 
instantaneous recording of data and its reduction, for imme- 
diate analysis, saves an immeasurable amount of time and 
reduces the chance of human error. The present experimental 
apparatus has been used over a range of Reynolds numbers 
throughout the laminar and transition regimes, into the tur- 
bulent regime. Results have compared well with previous data 
and correlations. A new flow channel of similar design should 
be constructed prior to further testing, due to the age and 
condition of the present equipment. A new set of calibrated 
thermocouples and possibly the introduction of pressure trans- 
ducers, instead of manometers, that could be connected to the 
data acquisition system, would provide more consistent, and 
as a result, more accurate data for flow rate measurements. 
Also, since the present apparatus has been modified several 
times in the past, resulting in leaks that are difficult to 
locate and stop, a new flow channel would help to eliminate 
some of the inaccuracies to temperature and flow measurements 
introduced by undiscovered leaks. It is also highly recom- 
mended that a method of safeguarding each piece of equipment, 
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so that in the event of a power failure, when power is re- 
stored, each item would require the operator to manually 
restart it . 
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APPENDIX A: EXPERIMENTAL UNCERTAINTY 



The uncertainties for the major variable in the experiments, 
were calculated in accordance with the method described by S. 
Kline and F. McClintoch [Ref. 29]. The estimates of the un- 
certainty in the measured quantities were made conservatively. 

As a result, there is considerable confidence in the uncer- 
tainties as calculated. 

The following equations were used to calculate the 
uncertainties : 



2 2 2 
(1) dA = A [YKA (^0 + (“) + (^) ) 

. , ,1/2 1/2 
dP a ; r 2 ,.p 2 

* (1/2 P air AP ((-Jii) ♦ ) ] 

air 



.2 d C - 2 d(T -T . ) 

r\ -r* o n -v nn t i -n ^ 



1/2 



(2) dQ • = Q • C(^) + ( - L a i- r ) + (- 

^ x air x air m v C • 

pair 



out in' 

T -T • 
out in 



0 ) 



(3) dh = h (( 



dQ air, 2 A . f dAT, 2 



1/2 



,t- 2 dD, , 2 dK . 2 

(4) dm = Nu ((— ) + c— + (-/^) ) 

TT D hd K air 



1/2 



, . 2 dD, , 2 dy . 2 dA 2 

(5) dRe hd - Re hd ( (£) + (-^) ♦ - C^) ) 

hd air c 



1/2 
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The values obtained for the uncertainties 



for the values 



in the curved 


section test run at a Reynolds number of 



17500 are: 



Quantity 


Uncertainty 


A 


.0005 


a pl 


.0026 


A 

c 


.0098 


c p 


.0034 


D hd 


.0087 


K 


.0412 


K 


. 0050 


K air 


.0012 


Nu 


.0421 


m 


.0223 


^air 


.0382 


Re hd 


.0260 


T blk 


.0092 


T. 

m 


.0046 


T 

1 out 


.0049 


T 

wo 


.0063 


T - T . 

out m 


.0136 


Y 


. 0051 


AP 


. 0222 


AT 


.0032 


p air 


.0027 


p air 


.0003 
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The major source of uncertainty in the average Nusselt number 
is the uncertainty in the pressure readings from the mano- 
meter which affects the mass flow rate calculation. In 
order to obtain a more accurate pressure measurement, it 
has been recommended that pressure transducers or inclined 
manometers be used instead of vertical manometers in future 
studies . 
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APPENDIX B: SAMPLE CALCULATIONS 



Figure 18 below, shows the major heat transfer components 
for each of the test sections. The sample calculations that 
follow, demonstrate the methods used by the computer to calcu- 
late these components as well as the Reynold's number, average 
heat transfer coefficient, and average Nusselt number for each 
set of data. The sample calculations are for the curved sec- 
tion, but those for the straight section are similar. Reynolds 
number, Nusselt number, and Dean number are based on hydraulic 
diameter . 




Figure 18. Energy balance in straight section. 
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A. SAMPLE CALCULATION DATA 



V H 


=: 


55. 2325 V 


V R 


= 


4.4075 V 


r pr 


= 


2 .0262 9 , 


p 

atm 


= 


30.45 in. Hg 


AP 


= 


22.5 in. H 9 0 


p i 


= 


13.05 in. Hy 


T in 


= 


25.41 C 


T 

out 


= 


28.08 C 


T 

wo 


= 


44.50 C 


T • 

Wl 


= 


25.99 C 


^insl 


= 


38.66 C 


^ins2 


= 


32.52 C 


T orif 


= 


28.32 C 


K air 


= 


26.52 x 10‘ 3 


'"pair 


= 


1006 J/Kg K 


^air 


= 


18.63 x 10' 6 


K ins 


= 


3.6495 x 10 


AX- 

ins 


= 


0.00635 m 


£ wo 


= 


0.70 


£ wi 


= 


0.40 


0 


= 


5.669 x 10" 


3 


= 


.5325 


Y 


= 


1.402 


gc 


= 


1 Kg m/NS 2 


R 


= 


286.8 Nm/Kg 



W/mK 

Kg/ms 
2 W/mK 



W/m 2 K 4 
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1.0 



wo - wi 



PL 



A 



A . 
pipe 

D c 

D pipe 

D • r 
onf 

p 

wet 



.07188 m 2 
.0016 m 2 
.002027 m 2 
.00635 m 
.0508 m 
.027051 m 
.5207 m 



B. TEMPERATURE CALCULATIONS 

1. Bulk Temperature 



T 



blk 



T . + T 

in out 

2 



25.41 + 28.08 
2 



26.75 C 



2. Mean Temperature Difference (AT) 



AT = T -Tv., = 44.50 - 26.75 = 17.75 C 
wo blk 

3. Temperature Difference in the Insulation (AT^ ng ) 



AT- 

ins 



^insl 



- T. 0 
ins 2 



= 38.66 



32.52 = 6.14 C 



C. POWER CALCULATIONS 

1. Power Supplied (Q ) 
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Q = _y — R = ( 55 /.. 2 2 2 .5H 4 / 4075 ) = 120.14 W 

P R PR (2.0262) 



2. Heat Lost Through the Outer Plate (Q^ 0 ) 



l0 " 



Qio = 

= 2.54 W 



AT- 

ins 



6.14 



(0.00635)/(3. 6495x10'-) ( 



Heat Radiated (Q ) 
a. Radiation Resistance (R^) 



R = l_fwo. I 

R A PL e „o 



1 - £ 



W1 



A p 

PL wo-wi PL wi 



1 r 1 , 1 11 - 1 r 1 „ 1 -n 

*pl e wi 1 ' Tior7my [ irT + TrT ^ 



b. Heat Radiated (Q ) 



Qr “ 



^ T wo 4 - T wi^ 



R 



T [= ] K 



R 



= (5. 669x1 0" 8 ) [(44.50 + 273) 4 - (25 .99 + 273) 4 ] 

40.74 



= 3.02 W 



D. MASS FLOW RATE CALCULATIONS 
1. Pressure Conversions 

a. P atm = ^" 48 x ^374.1 = 102741.35 N/m^ 



.07188) 



= 40.74 
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b. AP = 22.5 in H 2 0 x 248.64 = 5594.4 N/m“ 

c. P n = (13.05 in H ? 0 x 248.64) + 102741.35 = 105986.10 

1 7 

N/m 



2. Density of Air (p_. ) 

all 



atm 



air 



R T 



or if 



(102741.35) 

(2'8'6'.''8')'(Z873TT27Tr 



= 1.1889 Kg/m^ 



3. Expansion Factor (Y) 

Y = 1- [0.333 + 1. 145 (B 2 + 0.7B 5 + 12B 13 )] |L- 

Yi l 

= l-[0.333+l.l45((.5325) 2 +0.7(.5325) 5 +12(.5525) 10 ] 

5594.4 

x (1.402) (105986.10) 

= .9738 

4. Area of Orifice (A) 

A - li!W' . »C. 0. 2 . 7 0 5 J 2 . , 0005747 m 2 
4 4 

5. Mass Flow Rate (m) 

m = YKA /2gcP . XV = (. 9738) K( . 0005747) 

3. i r 

X /2(1) (1 . 1889) (5594 . 4) 
m = .06455 K 
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I terat ing : 

Assume a Reynolds number ^ = 55000 
' pipe 

Obtain a value for K, the flow coefficient, from 
reference 30. K = .6332 
Solve for m. m = .0409 Kg/s 

Solve for new Re - = ™ ^pipe = (. 0409) (. 0508) 

P P A p'ipe u 'air 1-002027) (18.63x10-6 



Re - = 54988 
pipe 

Check convergence and repeat 
(Convergence if difference 1 
m = .0409 Kg/s 



process 
ess than 



if necessary. 
. 001 ) 



E. REYNOLDS NUMBER CALCULATIONS 



m D 



1. Re 



'pipe 



P x P e 



A • y - 
pipe air 



( .0409) ( .0508) 
CTW2T3T77 (l8.63xlU- 6 



54988 



0 n m D c ( .0409) ( .00635) 

• Ke d A c y— ( .0016) ( 18.6 5x1 U * 6 



8713 



m D 



3 . Re 



hd 



hd 



A y - 
c air 



a . 




4 * A c 
^wet 



(4) ( .0016) 

rrnn 



.01229 m 



m D 



Re 



hd 



hd 



A c ^air 



( .0409) ( .01229) 
(.0016) (18.63x10-6 



16863 



F. HEAT CONVECTED TO AIR CALCULATION 

1. Q air = ^Cp^out^in^ = (.0409) (1006) (28. 08-25. 41)-109.86W 
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G. AVERAGE HEAT TRANSFER COEFFICIENT CALCULATION 



h = 



'air 



109 . 86 



A pl AT 



(.07188) (17.75) 



= 86.11 W/m C 



H. AVERAGE NUSSELT NUMBER CALCULATION 

m = K D hd = (86.11) (.01229) = 39>go 

K • 26.52 x 10'- 5 

air 



I. DEAN NUMBER CALCULATION 



De 



Re 



hd 




16863 



01229 
. 505 



33 85 
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APPENDIX C: CORRELATIONS 



The correlations obtained for the present study were 

obtained using the method of least squares for a first degree 

polynomial, as outlined by C. F. Gerald [Ref. 29]. In this 

method the values of Reynolds number (x) and Nusselt number 

(Y) were first converted to their natural logarithmic value. 

2 

Next the summations of x^, x^“ , Y^, and x^Y^ were calculated 
The values of these quantities were then placed in a matrix: 



[ N Ex i 2 ] [ 1 ° } = [ S Yi ] 

Ex^x^ 1 Ex i Y • 



and the resulting simultaneous equations 



a N + a-. Ex • = E Y,- 

O 11 1 



a o zx i * a l Zx i = Zx i Y i 



were solved for a^ and a n . 

o 1 

The equation of the line then became: 

Ln Nu = a 0 + a-^Ln Re 
or 

m = e a ° Re ai 

The standard error for the present study correlation was 
calculated by the following equation: 

2 E(Y i - 
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where: is the actual value of the average Nusselt number 

obtained experimentally. 

y^ is the value of the average Nusselt number calculated 
from the correlation equation. 

N is the number of data points used in the correlation, 
n is the degree of the polynomial used. 

For the present study the following correlation was obtained 
for the straight section: 

Nu = .063 Re hd °’ 65 

with a standard error of: 
a 2 = 0.1988 

The correlation obtained for the curved section was: 

Nu = .040 Re hd °’ 72 

and its standard error was determined to be: 
a 2 = 0.1050 
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